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A model-based description of the scaling and radial location of turbulent fluctuations 
in turbulent pipe flow is presented and used to illuminate the scaling behaviour of the 
very large scale motions. The model is derived by treating the nonlinearity in the pertur- 
bation equation (involving the Reynolds stress) as an unknown forcing, yielding a linear 
relationship between the velocity field response and this nonlinearity. We do not assume 
small perturbations. We examine propagating modes, permitting comparison of our re- 
sults to experimental data, and identify the steady component of the velocity field that 
varies only in the wall-normal direction as the turbulent mean profile. The "optimal" 
forcing shape, that gives the largest velocity response, is assumed to lead to modes that 
will be dominant and hence observed in turbulent pipe flow. 

An investigation of the most amplified velocity response at a given wavenumber- 
frequency combination reveals critical layer-like behaviour reminiscent of the neutrally 
stable solutions of the Orr-Sommerfeld equation in linearly unstable flow. Two distinct 
regions in the flow where the influence of viscosity becomes important can be identified, 
namely a wall layer that scales with and a critical layer, where the propagation 

velocity is equal to the local mean velocity, that scales with R^^/^ in pipe flow. This 
framework appears to be consistent with several scaling results in wall turbulence and 
reveals a mechanism by which the effects of viscosity can extend well beyond the imme- 
diate vicinity of the wall. 

The model reproduces inner scaling of the small scales near the wall and an approach 
to outer scaling in the flow interior. The appropriate scaling velocity for the very large 
scale motions is predicted to be the centreline velocity, in agreement with experimental 
results. Lastly, we interpret the wall modes as the motion required to meet the wall 
boundary condition, identifying the interaction between the critical and wall modes as 
a potential origin for an interaction between the large and small scales that has been 
observed in recent literature as an amplitude modulation of the near-wall turbulence by 
the very large scales. 



1. Introduction 

The accurate description of statistical scaling and instantaneous structural coherence 
of turbulent fluctuations, and their relationship to the mean flow, arc amongst the im- 
portant unsolved problems in physics. Even at the simplest level, we still lack a complete 
explanation of the development of the mean velocity profile in canonical flows as the 
Reynolds number increases. Such an understanding would be an important step towards 
accurate prediction of skin friction in complex aeronautical and industrial flows, and 
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would underpin any effort at a turbulence control scheme. A brief summary of three 
distinct approaches to these problems is given here. 

The research literature is replete with statistical descriptions of the wall-normal dis- 
tributions of the components of the Reynolds stress tensor, with an obvious bias to the 
streamwisc normal stress, which is most easily measured. New understanding of turbu- 
lence has emerged in the past decade. This has included differences in some characteris- 
tics between the canonical cases of pipe, channel and boundary layer flow (Monty et al. 
2009), in contrast to the long-standing hypothesis of universality of near- wall scaling. 
Relevant to the present work, several studies have revealed highly energetic structures 
with streamwise wavelength of order ten times the outer lengthscale, deemed variously 
Very Large Scale Motions (VLSMs) or superstructures. The VLSM phenomenon, dis- 
cussed at length later in this work, suggests a more complicated nature of the scaling 
of the turbulence than an inner/outer/overlap structure proposed for the mean velocity. 
This has been confirmed to the resolution of experimental measurements, for the mean 
velocity. Specifically, it confirms the influence on the inner, near-wall flow, of turbulent 
activity that scales on outer variables. 

The advent of Particle Image Velocimetry (PIV) techniques and advances in simulation, 
alongside more traditional visualisation techniques, have illuminated the development 
and grouping of organised, coherent structures. These include the autonomous near-wall 
cycle and an eddy type that is statistically well described by the hairpin vortex paradigm. 

In parallel to these experimental and computational approaches, considerable progress 
has been made in understanding the amplification properties of the Navier-Stokes oper- 
ator, and in particular the linearised operator, in laminar flows. More recent work has 
extended some of these techniques to the turbulent case, with limited success. 

In what follows, we summarise some key concepts and questions arising from these 
three distinct approaches which are pertinent to the current work. 

1.1. The challenge to classical scaling: the influence of VLSMs on the near-wall region 

Classical scaling ideas involve inner and outer layers, where the appropriate spatial scales 
are respectively the viscous unit and the outer lengthscale, and the appropriate velocity 
scale is the friction velocity. In the case of the outer layer this velocity is impressed by the 
boundary condition at the wall. For sufficiently high Reynolds number, there may also 
be an overlap layer in which both scalings hold and therefore the important lengthscale 
must be the local distance from the wall. This scaling appears to work well for the mean 
velocity despite recent challenges (Marusic et al. 2009). However it is unable to collapse 
the turbulent fluctuations in the inner and overlap regions, an apparent reflection of the 
influence of outer scales on the inner region, termed "inner-outer interaction" as explored 
by Bandyopadhyay & Hussain (1984), that is reflected in the spectral energy distribution. 

The recent focus on VLSMs in canonical turbulent wall flows has given some insight 
into the source of this interaction. Coherence across the wall layer at long streamwise 
lengthscales has been known since the observations of Kovasznay et al. (1970). The work 
of Kim & Lim (1993) and Morrison et al. (2004) identified the energetic importance 
of very large scale motions in the streamwise spectra in turbulent pipe fiow and the 
subsequent work by Adrian and co-workers (Guala et al. 2006; Balakumar & Adrian 2007) 
has indicated that large and very large scale features must be considered to be "active" 
in the sense that they carry significant shear stress. This contrasts with Townscnd's 
geometrical arguments that only the small scales display sufficient coherence in the wall- 
normal and streamwise fluctuations to extract signiflcant energy from the mean flow, via 
the product of the Reynolds shear stress and the mean velocity gradient. 

Simulating and observing these high aspect ratio VLSMs places heavy demands on 
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existing computational and experimental techniques, a constraint that has obstructed 
progress in our understanding of their origin and development in turbulent flows. The 
statistical imprint in the streamwise and wall-normal directions is clear from hot-wire 
measurements. Hutchins & Marusic (20076) showed for sufhcicntly high Reynolds num- 
bers that this imprint reaches from a peak energy located somewhere in the overlap re- 
gion down to the wall, even having a footprint on the wall shear stress (Marusic & Heuer 
2007): Hutchins & Marusic (2007a) and Monty et al. (2007) used Taylor's hypothesis 
to reconstruct the long streamwise-spanwise coherence in the streamwise velocity from 
arrays of hot-wires in turbulent boundary layers and pipe flow, respectively. Their use of 
Taylor's hypothesis highlights a major observational difficulty: at these large scales, Tay- 
lor's hypothesis may not apply, since the necessary arguments based on the ratio of eddy 
turnover timescale to convective timescale no longer hold. In addition, the wall-normal 
extent of the VLSMs means that if the VLSM convects with the local mean velocity 
somewhere in the overlap layer, Taylor's hypothesis with a scale-independent convective 
velocity must be in error closer to the wall, and likely also further from the wall to a 
much lesser degree. 

Having listed some of the outstanding questions concerning the VLSMs, we now provide 
a brief summary of what is known about their extent and influence on the near-wall 
region. Following the first observation by Hutchins & Marusic that the peak streamwise 
VLSM energy appeared to occur at a constant fraction of the boundary layer thickness, 
y/5 ~ 0.05, the location of the peak has been shown by McKeon (2008) and Mathis et al. 
(2009) to have a weak Reynolds-number dependence when a sufficiently large range of 
Reynolds number considered. The latter authors compared the peak's location with the 
centre of a mean velocity overlap region that has an inner limit that is either fixed in 
inner units or Reynolds number dependent, giving rise to 5'^^^'^ and S~^^^^ dependence of 
the peak location, respectively, where S'^ = Sut/v, v is the kinematic viscosity and the 
friction velocity Ur — \/t j p. However the agreement depends crucially on single data 
points obtained at very high Reynolds number in the near-neutral atmospheric surface 
layer, and these vary between studies (Mathis et al. 2009; Metzger et al. 2007; Guala 
et al. 2009 &). 

Recent studies show that the streamwise extent of the dominant large scale motion 
is proportional to the outer flow lengthscale and is closely given by Xx,vlsm lOR 
in internal flows (Kim & Adrian 1999; Monty et al. 2009) and Xx, superstructure ^ 6(5 in 
boundary layers (Hutchins & Marusic 2007&), while the Direct Numerical Simulation 
(DNS) of channel flow of Jimenez et al. (2004) and the study of Monty et al. indicate 
the additional importance of a slightly smaller wavelength, Xx,lsm ^ 3/i. Here R, 6 and 
h are the pipe radius, boundary layer thickness and channel half-height, respectively. 
Monty et al. (2009) recently performed a rigorous comparison of the streamwise velocity 
spectra in the different flows at the same Reynolds number and elucidated the respective 
importance of these three wavelengths in different regions of the flow. The spanwise 
extent in each case appears to be of the order of one outer lengthscale (Hutchins & 
Marusic 2007a; Monty et al. 2007), giving the largest scale structure an approximate 
axially elongated aspect ratio of 10 : 1 : 1. 

Subsequently Mathis et al. (2009) expounded on earlier observations that the very large 
scales apply an amplitude modulation to the small scale turbulence near the wall. Using a 
Hilbert transform technique, they were able to quantify the interaction and its Reynolds 
number dependence, as well as demonstrate a change in the sign of the modulation that 
corresponded with the location of the VLSM energy peak. Guala et al. (2009a), Guala 
et al. (20096) and Chung & McKeon (2009) have also proposed methodologies to describe 
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this effect, with the latter identifying that the modulation can also be described in terms 
of the spatial phase relationship between large and small scale turbulent activity. 

It is clear, then, that the very large scale structure reflects an aspect of boundary layer 
dynamics that has hitherto been poorly understood, but has importance for scaling of 
local and global turbulence properties, as well as for future flow control schemes. 

1.2. Alternative scaling approaches from theory and observation: Critical Layers 

Several theories have been proposed to account for the missing physics in the classical 
scaling, including the mesolayer of Long & Chen (1981), the focus on the location of the 
Reynolds stress peak by Sreenivasan and co-workers (Sreenivasan & Sahay 1997) and 
the hierarchical structure associated with the mean momentum balance of Klewicki and 
co-workers (Klewicki et al. 2007). 

Sreenivasan (1988) formulated an inviscid structural model of the turbulent boundary 
layer consisting of two vortex sheets of opposite sign symmetrically located about the 
hypothetical wall location. By analogy with the critical layer in transitional boundary 
layers, he showed that the wall-normal location of the peak in Reynolds shear stress 
should have a Reynolds number dependence given by y'^^ ^ _R+^/^ with a proportionality 
constant of two determined from experimental data. Here y+ = yur/v is the wall-normal 
distance y non-dimensionalised with the viscous scaling length, defined as the kinematic 
viscosity v divided by the friction velocity. The wall shear stress is denoted Tw and p 
represents the density. Experimental evidence from boundary layers, pipe and channel 
flows supported this critical layer interpretation, the resulting streamwise and spanwise 
wavelengths for maximum amplification, the y'^^, ^ scaling (Sreenivasan & Sahay 

1997) and the hypothesis that the mean velocity at the peak in the Reynolds shear stress 
corresponds to a constant fraction of the freestream velocity. In addition, Sreenivasan 
& Bcrshadskii (2006) were able to extend the study of the region in the vicinity of the 
Reynolds stress peak by performing a logarithmic expansion and predicting the resultant 
form for the mean velocity profile. 

The reader should note that an analogy with critical layer theory in laminar flow 
was also made by Sirovich et al. (1990), who found that the most energetic propagating 
eigenfunctions in a Proper Orthogonal Dccompositionf (POD) of turbulent channel flow 
had their principal support in the region of the peak in Reynolds shear stress. Duggleby 
et al. (2007) have also given some insight into the most energetic POD modes in low 
Reynolds number turbulent pipe flow, illuminating a distinction between propagating 
and non-propagating, wall, lift, asymmetry and ring modes defined by the relative mag- 
nitude of the streamwise and azimuthal wavenumbers. A truncated POD representation 
of turbulent pipe flow was also used by Aubry et al. (1988) to develop a model for the 
dynamical behaviour of the streamwise roll modes and establish a connection between 
near-wall turbulent flow and the dynamics of a chaotic system. 

1.3. Linear and nonlinear tools for analysis of the Navier-Stokes equations 

Other researchers, essentially forming a different community, have exploited tools from 
linear systems and dynamical systems theory, with a notable focus on the large algebraic 
energy growth that is possible due to the non- normality of a stable Navier-Stokes operator 
that has been linearised about a base flow. In this picture of transition, for sufficiently 
large perturbations, such growth is understood to then bring the locally stable system out 
of its basin of attraction, inducing nonlinear behaviour that is associated with turbulence. 
An operator A is non-normal if AA* ^ A* A where A* denotes the adjoint of A. The 

f also known as Karhunen-Loeve analysis 
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adjoint is defined with respect to an inner product, so we see that non-normahty is only 
defined with respect to a particular inner product. Typically, for this type of study, the 
C2 or perturbation energy norm is of interest. Flows that are linearised about a steady 
flow solution with shearing can yield highly non-normal operators. In the past fifteen 
years, significant progress has been made in understanding system non-normality as a 
mechanism for energy amplification in shear fiows (Butler & Farrell 19926; Trefethen 
et al. 1993; Farrell & loannou 1993) leading to nonlinear breakdown in both linearly 
stable and unstable fiows (Jovanovic & Bamieh 2004). 

There has been less investigation of the non-normal growth mechanisms in turbulent 
fiow, with the notable exceptions of the early attempt by Butler & Farrell (19926) to 
predict the spacing of near-wall streaks in turbulent flow and the more recent studies 
of del Alamo & Jimenez (2006), Cossu et al. (2009) and Wilhs et al. (2009). A notable 
and somewhat limiting issue in the treatment of turbulent fiow is the modelling of the 
interaction of the amplified disturbances with the "background" turbulence. A solution 
introduced by Reynolds & Hussain (1972) has been to use the eddy viscosity formulation 
of Cess (1958), but this relies on an a priori knowledge of the spatially-averaged, wall- 
normal variation of the mean Reynolds stress integrated across contributions from various 
Reynolds numbers. Other attempts have been made to use linear analysis to explain the 
dominant features of turbulent flow in terms of optimal transient modes in the initial value 
problem (Butler & Farrell 19926; del Alamo & Jimenez 2006; Cossu et al. 2009), response 
to stochastic forcing (Farrell & loannou 1998; Bamieh & Dahleh 2001) and system norm 
analysis (Jovanovic & Bamieh 2005; Meseguer & Trefethen 2003). In recent work, Willis 
et al. (2009) have investigated the maximal response to harmonic forcing in pipe flow. 
Perhaps most importantly, it has been shown that both linear non-normality (Hcnningson 
& Rcddy 1994) and the terms that are linear in the turbulent fluctuation (Kim & Lim 
1993) are required to sustain turbulence in infinite or periodic wall-bounded fiows. 

del Alamo & Jimenez (2006) made a direct comparison between a transient growth 
analysis of channel flow and earlier DNS results, and showed that while the analysis could 
predict two spanwise wavelengths that would experience large transient energy growth 
that were in good agreement with the computational and experimental observations, 
the corresponding predicted streamwise wavelengths were significantly too high. Cossu 
et al. (2009) performed a similar analysis in a turbulent boundary layer, modelling it as 
a parallel flow. The realisation that turbulent wall flows can be linearly stable followed 
the work of Reynolds & Tiederman (1967), who demonstrated this for turbulent channel 
flow. 

In terms of the physical mechanisms and structure behind the energetic near- wall cycle, 
Schoppa & Hussain (2002) used linear perturbation methods to propose a streak transient 
growth mechanism capable of reproducing with good fldelity structures observed in a 
low Reynolds number DNS of a channel. Waleffe (1997, 2001, 2003) followed a different 
approach, developing exact solutions of the Navier-Stokes equations that give rise to 
unstable coherent structures which they propose form the foundations of transitional fiow 
and near-wall turbulence. Subsequently there has been much interest in the importance 
and observability of travelling wave solutions of the Navier-Stokes equations, e.g. Wedin 
& Kerswell (2004) and Viswanath (2009) in pipe flow (which, of course, is also linearly 
stable). The ongoing work of Gayme et al. (2009) with a forced, two-dimensional, three 
velocity component model provides another attempt to predict the form of the turbulent 
mean flow. 
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1.4. Selection of pipe flow for further study 



Of the canonical turbulent flows, pipe flow has received considerable attention since 
Osborne Reynolds' seminal work identifying the role of the Reynolds number in the flow 
behaviour. Pipe flow is also important to many obvious industrial applications. From an 
experimental point of view, providing that the development length is sufficiently long 
for the flow to be considered to be fully-developed, a pipe constitutes a well-defined, 
simple to generate geometry. As such, there is a wealth of experimental pipe fiow data 
available for comparison, with the results from the Princeton/ONR Superpipe providing 
detailed information on the mean velocity (Zagarola & Smits 1998; McKeon et al. 2004), 
streamwise (Morrison et al. 2004) and wall-normal (Zhao & Smits 2007) fiuctuations, and 
azimuthal correlations (Bailey et al. 2008; Bailey & Smits 2009) across three decades 
in Reynolds number. The results of Monty et al. (2007) at an intermediate Reynolds 
number have already been described in the context of the VLSMs. In addition, a recent 
DNS study by Wu & Moin (2008) explored the properties of the spatial velocity field 
at relatively low turbulent Reynolds number, while the POD analysis of Duggleby et al. 
(2007) categorised energetic propagating and non-propagating modes a posteriori, from 
full field information from a DNS at i?+ = 180. Lastly, the geometry of the pipe has the 
useful property of imposing a restriction on the azimuthal wavenumber to integer values 
only, which simplifies our analysis in the subsequent sections. 



In this work we propose a simple analysis of the Navier-Stokes equations for incom- 
pressible, fully-developed turbulent pipe fiow. The resulting model describes the spatial 
distribution of turbulent energy in the three-dimensional propagating velocity modes 
that are most responsive to harmonic forcing. Wc offer the model as a first step towards 
bridging the gap between statistical and structural interpretations of wall turbulence, 
since it provides qualitative information on both the temporal and spatial distributions 
of velocity associated with each mode. We suggest that such a reconciliation of the differ- 
ing observations of the same system would provide an important advance in the field. In 
what follows, we use the model to investigate an issue of intense current interest in the 
boundary layer community; the radial extent, characteristics and scaling of very large 
scale motions. 



2. A simple model for the spatio-temporal distribution of turbulent 
energy in pipe flow 

In following analysis, we develop a formulation of the Navier-Stokes equations designed 
to examine the receptivity of turbulent pipe fiow to forcing. In this way, we develop 
a framework that permits investigation of the form and likely magnitude of spatially- 
and temporally-harmonic, propagating finite-amplitude fluctuations about the turbulent 
mean profile in pipe fiow. 

2.1. Pipe flow equations and non-dimensionalisation 

The non-dimensional Navier-Stokes equations for fully-developed, incompressible pipe 
flow with constant viscosity arc given by 



1.5. Contribution 



dtu 



Vp - u • Vu + — V^u 



Re 



(2.1) 
(2.2) 



V-u=0 



A critical layer model for turbulent pipe flow 



7 




Figure 1. Schematic of pipe geometry and nomenclature. 



We follow the eonvention of Meseguer & Trefethen (2003), where the equations of motion 
are non-dimensionalised with respect to the pipe diameter and twice the bulk, volume- 
averaged velocity, Ubuik (which in their study is equal to the laminar flow centreline 
velocity). Thus the Reynolds number in Equation 2.2 can be defined as 

„ UbulkD 

Re = . 

V 

Here D is the pipe diameter and v is the kinematic viscosity. We retain the boundary 
layer terminology in fixing y = 1 — r, and w, v(= — w') and w as corresponding to the 
streamwise, wall- normal and azimuthal velocities such that u = {v',w,u), as shown in 
figure 1. 

2.2. Model development 

We introduce a projection onto a divergence free basis: {Cm(^)} in the radial direction, 
Fourier modes in the homogeneous spatial directions and the Laplace transform in time. 
Implicitly, we are considering a pipe that is infinitely long or periodic in the axial direc- 
tion. Assuming fully developed flow allows us to express the velocity field as the sum of 
harmonic modes. Then, 

u{r,x,9,t)^—Y, / c™fc„^e™(Oe^'"+"''+'*dfcds (2.3) 



27ri ^ 



?'00 pOG 

Cm k n LO ^rn k n uj 

(r)dfcds (2.4) 



%oo J —oo 



with s ~ iuj, so that only harmonic forcing and response is considered. 

The wavenumbers (fc, n, lu) are non-dimensional such that k = k' R, etc. The integration 
path for the Laplace transform is over the closed right half plane, which is analytic in the 
case of pipe flow. In the inviscid limit special treatment for singularities on the imaginary 
axis would be required, but we do not consider this limit. 

We work with the spatial £2 inner product throughout. 



{a,b) ~ / a{x)b{x)dx. 
Jxefi 



(2.5) 
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The basis functions are required to have the special properties 

{^aAb)=5at, (2.6) 

V ■ ^ir)^knu, = 0. (2.7) 

This is done to ehminate the pressure term. For notational convenience we make the 
definitions 

Ufe„. = (u,e'("*+''^-+"'')) (2.8) 

u = u - uooo (2.9) 

fknu = (-U ■ Va, e'(-*+fe-+»«)) . (2.f 0) 



This finaUy yields equations for the fluctuations that are linear in Uknuj, and a base 
flow equation, 

iUJUknui = ^knuUknui + ^knui, V(fc, 71, w) ^ (0, 0, 0), (2.ff) 

= fooo - Uooo ■ Vuooo + ^V^uooo. (2.f2) 

Re 

The unknown constant forcing fooo describes the maintenance of Uooo via the radial 
derivative of the Reynolds stresses, generated from interaction with the other modes. 
We can identify Uooo with the turbulent mean velocity profile. Similarly, iknuj describes 
the excitation of Uknuj by the triadic interaction with other wavenumbers. We cannot 
solve these equations without additional information because Cknu incorporates Uooo and 
fooo incorporates u. This is essentially an appearance of the closure problem. However, 
the reader will notice that the perturbation equation (2.11) is a linear system with an 
unknown forcing iknu- This fact is central to our treatment. 

We avoid the closure problem for the base flow equation simply by knowing a priori the 
mean profile from experimental data. This allows calculation of Cknuj and precludes the 
need for an eddy viscosity formulation normally required, e.g. by (Reynolds & Hussain 
1972)del Alamo & Jimenez (2006), etc. 

While fooo can be simply calculated from the mean profile, we do not know the Reynolds 
stress at any other wavenumber combination, tknio- Our approach is not to make as- 
sumptions about this forcing, but to simply examine the response of (2.11) at individual 
{k,n,ui) triplets over the set of all possible harmonic forcings. 

The linear operator has the explicit form for pipe flow 

/ Re-\d^ +r-^dr -n'^r-'^ - k"^ -r-'^) -ikuQoa - 2inr-'^ 

jOkriLj = 2mr^^ Re^^{d^ + r^^dr — n'^r^'^ — k"^ — r^'^) - ikuooo 

y — 9,.uooo Re~^{d^ + r~^dr — n'^r^^ — fc^) — ifcwooo 

(2.13) 

where, as before, the states in u are the radial, azimuthal and axial velocities respectively. 
We make the following observations: 

(a) Only integer n are permissible; 

(&) We expect downstream travelling waves such that k and u arc of opposite sign: 
henceforth k > and w — 5- — w; 

While Equation 2.11 is linear in Uknuj, and Cknuj is identical to the operator ob- 
tained by linearising around the turbulent mean velocity profile, no linearisation has 
been performed. Nonlinear effects at other wavenumber-frequency combinations are re- 
tained through the action of the forcing fknui- In physical terms, this forcing can be 
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considered to stimulate fluctuations that may lead to a net energy gain because of the 
characteristics of Cknu ■ 

The model fully describes the energetic interaction between the base flow and fluc- 
tuation Uknui, given a mean profile. Additionally, fkuLu acts perpendicular to Uknu, that 
is, {fknui , Uknuj) = 0. This Can be interpreted as ffc™^ being conservative with respect to 
fluctuation energy, and so responsible for the transfer of energy in spectral space but not 
directly responsible for the extraction of energy from the base flow. A similar formula- 
tion to the current one, using this fact to derive a globally laminarising control law, was 
described in Sharma et al. (2006). 

Next, we proceed to analyse the response of CkuLo at a particular wavenumber combi- 
nation, subjected to the harmonic forcing iknu- 

2.3. Resolvent norms and model formulation 
Equation 2.11 can be rearranged as 

Ufcno; = (jw/ - CknuiY^hnu,- (2.14) 

The operator [iwl — Cknuj)^^ is called the resolvent and is the focus of our analysis. It 
provides a measure of the turbulent energy response that is possible for a given forcing. 
Some interpretation of the resolvent is given in Appendix A. 

For pipe flow, using (2.13) and for (fc, n, ur) ^ (0, 0, 0), the resolvent can be written as: 

ikU-iuj B 

{iujI-Cknu,y = \ -B - ARe-^ + ikU - iuj 

drU - ARer^ + ikU - iuj 

(2.15) 

with A = {d^ + \dr — n^r ^ ~ k'^ — r B = 2inr ^ and U — uqoo- Equation 2.15 
clearly highlights that the operator is not self-adjoint in the presence of drU. The shear 
is a source of non-normality under the energy norm, coupling the radial and axial velocity 
components. The A/ Re term will always be small for Reynolds number large enough for 
turbulent flow. One might expect large values of the resolvent norm under any of the 
following conditions: 

(a) in regions of high shear, where the drU is large; 

(6) at critical layers where uj/k = U{y), so the component of the normal speed of 
propagation of the wave in the streamwise direction is equal to the local mean velocity; 
(c) for stationary modes with fc = cj = 0. 

2.4. Most amplified modes 

We seek a decomposition of the resolvent at a particular wavenumber pair and frequency 
which ranks the response to forcing in some sense. We take the Schmidt decomposition 
(called the singular value decomposition in the discrete case) of the resolvent, namely 

oo 

{iuil - Cknu)^^ = ^il^j{k,n,y,uj)aj{k,n,uj)(l)*{k,n,y,uj) (2-16) 

with an orthogonality condition 

(j)i{k,n,y,ijj)(j)m{k,n,y,uj)Ay = 5im (2-17) 

'4>i{k,n,y,Lo)'4)rn{k,n,y,uj)Ay 5im (2-18) 

V 
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and 

o-i > o-;+i > 0. 

The and "i/jj form the right and left Schmidt bases for the forcing and velocity 
fields and the real aj are the singular values. This decomposition exists if there are no 
eigenvalues of L with zero real part and is unique up to a pre-multiplying unitary complex 
factor on both bases corresponding to a phase shift and up to the ordering of the cTj's 
(Young 1988). 

This basis pair can then be used to decompose arbitrary forcing and the resulting 
velocity at a particular Fourier component 

(2.19) 
(2.20) 

The energy of the same Fourier component of the resulting disturbance velocity is 

oo 
1=1 

Clearly the forcing shape that gives the largest energy at a particular frequency and 
wavenumber is given with aiknu} =0, 1^1. This approach permits the investigation of the 
dependence of maximum energy amplification on the form of the forcing in the wavenum- 
ber and frequency domain. The singular value decomposition for a given wavenumber pair 
and frequency corresponds to full volume, three component forcing and response modes 
ranked by the receptivity of the linear Navier-Stokes operator. The velocity response 
must have the same k and n but not necessarily the same y distribution (spatial phase 
variation in y) as the forcing. 

By Parseval's theorem, the energy integrated over frequency and wavenumber is equal 
to that integrated over the temporal and spatial domains (the spectral and physical 
spaces are isomorphic) . As such, the C2 norm of the resolvent is its leading singular value, 
(7i. This means that the normalised harmonic forcing that gives the largest disturbance 
energy in the C2{^ x [0, 00)) sense is f = (f>i, with a 'gain' of cti. The next largest 
arises from f = (t>2 and so on, at a particular wavenumber pair and frequency. The 
corresponding flow response modes are given by the related w = i/'jf , "02 1 etc. For aj near 
zero, the modes are not easily computed because they are effectively degenerate. However 
for the leading singular values the mode shapes are extremely robust to numerical error. 

This decomposition permits analysis of what we call the forcing and response modes 
associated with large responses of the flow (the "optimal response"). In this sense, this 
decomposition analyses the receptivity of the flow to forcing. 

In what follows, we focus our attention on the modes associated with the first singular 
value for a range of wavenumber-frequency combinations {k,n,uj) and show that they 
agree very well with experimental observations and classical scaling concepts. 

2.5. Computational approach 

The computational analysis of the linear operator, was performed using a modified 
version of the spectral code of Meseguer & Trefethen (2003). The code essentially provides 
the operator jCknui described above, evaluated at a number N of wall-normal grid points. 
For the particular problem under consideration here, the only modification to the linear 
operator used by Meseguer & Trefethen (2003) is the use of the turbulent mean velocity 
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profile instead of the steady, parabolic laminar base flow, as discussed above. We retain 
the same non-dimensionalisation, namely using the centreline velocity for a laminar flow 
with an equivalent mass flux, the steady laminar flow pressure gradient and the pipe 
radius, a. This formulation is equivalent to using the bulk-averaged velocity, Ubuik, and 
pipe diameter, D, as velocity and lengthscales, which is more natural to the turbulent 
problem. Of course, because the turbulent velocity profile is blunter than a laminar one, 
the constant mass-flux constraint means that the non-dimensional turbulent centreline 
velocity will always be less than one. 

The form of the turbulent mean profiles was determined directly from experimen- 
tal data obtained using Pitot probes in the Princeton/ONR Superpipe and reported 
by McKeon et al. (2004). This experimental data spans the Reynolds number range 
31 X 10^ <; i?e 35 X 10^ or 860 R+ < 535 x lO^, where Ur = V'Wp is the 
friction velocity, is the mean wall shear stress and p is the density. Interpolation be- 
tween Reynolds numbers was performed in a process equivalent to assuming a Rcynolds- 
number-independcnt form of the velocity profile. Issues of numerical stability limit our 
study to Reynolds numbers Re ^ 10 x 10^, well below the estimate of the Reynolds num- 
ber at which the Superpipe results may show some effect of wall roughness, Re ^ 24 x 10^ 
(Shockling et al. 2006). Note that this Reynolds number range also spans conditions that 
have been described as representative of high Reynolds number turbulence in pipe flow 
by McKeon & Morrison (2007) (i?"*" > 5 x lO'^) and in boundary layers by Hutchins & 
Marusic (20076) (6+ > 4 x 10^). 

The very sensitivity of the operator to perturbations that leads to large energy am- 
plification can also cause numerical stability issues associated with the spatial resolution 
employed in the pipe cross-section. Judicious choice of this resolution, A^, is required. 
This is briefly discussed in Appendix B. 

2.6. A model of spectral energy distribution 

While the concept of the Navier-Stokes equations as a linear system with a non-linear 
feedback forcing was been raised in a control theory context by Sharma et al. (2006), the 
emphasis here is more toward the receptivity of the flow, as explored for laminar pipe flow 
by Sharma & McKeon (2009). This approach is comparable to the "mother-daughter" 
scenario of Boberg & Brosa (1988), where a linear but non-normal process allows small 
disturbances to feed more energetic disturbances that can dissipate the gained energy. 
Non-linear effects then transfer some of this energy to the smaller initial disturbances. 
The understandings differ to the extent that the "mother-daughter" scenario considers 
the evolution of structural perturbations, and naturally leads to the transient growth 
problem of Butler & Farrell (1992&) and others. In contrast, the current work considers 
the gain response to harmonic forcing, naturally leading to a linear input-output analysis. 

The modes under investigation are propagating in the streamwise and spanwise direc- 
tions and distributed in the wall-normal direction. As such they are strongly analogous 
to the spectral decomposition of spatial and temporal velocity flelds from experiments 
and simulations. Their wave-like nature implies that the propagation (phase) velocity of 
each mode is given by = uj/K, where K = {k'^ + n'^y/^, with str eamwise (normal) 
component Ux = uj/k. In global terms, we expect that this analysis will give some in- 
sight into the spectral energy storage. Under forcing of the appropriate magnitude at all 
wavenumbers and frequencies, the correct Reynolds stress tensor would be calculated, 
obtaining both the true variation of turbulent energy production and dissipation, and 
closure of the feedback to the given turbulent mean velocity proflle. 

In this spirit, the current work explores the (fc, n, lo, Re) parameter space and validates 
the efficacy of our simple model, by demonstrating that it is capable of capturing classical 



12 



B. J. McKeon and A. S. Sharma 



and empirically-observed features of wall turbulence, such as inner and outer scaling, 
"attached" and "detached" motions and Reynolds number trends. 

3. Predictions of the model 

In this section, we describe the characteristics of the modes that are predicted by 
the preceding analysis. We consider only propagating, helical modes in the streamwise 
direction, with fc 7^ 0, n 7^ 0. 

We begin by presenting results for Reynolds numbers Re = 75 x lO'^ and Re = 410 x 10'^ 
= 1800 and 8500, respectively). The former is representative of the upper range 
accessible with current DNS techniques, while the latter Reynolds number is high enough 
to ensure characteristics of "high Reynolds number" turbulence (McKeon et al. 2004; 
McKeon & Morrison 2007; Hutchins & Marusic 20076). In addition, we know the velocity 
profiles and the streamwise and wall normal turbulence statistics and spectra from the 
Superpipe at, or close to, these conditions. This facilitates comparison with experimental 
results (McKcon et al. 2004; Morrison et al. 2004; Zhao & Smits 2007). From these data, 
we expect that the streamwise turbulent energy across the pipe cross-section lies within 
the range 5 x 10"^ < < 5 ^ 102. 

3.1. Perturbation mode shapes 

Figure 2 compares the velocity distributions for modes with the wavenumber pair (fc, n, lu) = 
(2, 2, 1) for the first three singular values, am with m = 1,2, 3. Clearly the higher order 
modes generate velocity distributions with more local maxima in the radial direction, 
apparently m maxima in the streamwise velocity for the modes shown. In addition, the 
singular value Bode plot of figure 3 shows that the first singular value is significantly 
larger than the other singular values. This suggests that the mode corresponding to ci 
is likely to dominate in observations in real pipe fiow. By comparison. Bailey & Smits 
(2009) have shown that more than 75% of the streamwise energy is contained in the first 
radial POD mode for Re — 150 x 10^. Thus we focus on the the first mode. 

The literature associated with the amplification properties of the linearised Navier- 
Stokes operator indicates that maximum amplification of energy is obtained for the 
streamwise constant modes. While we do not consider a fc of exactly zero in this study, 
we do observe similar trends as fc — 0. However we focus instead on the form of the ve- 
locity modes with streamwise wavenumber in the range corresponding to that observed 
in experiments. 

Figures 4-6 show the distribution of turbulent energy in the pipe for the first mode with 
wavenumber pair (fc, n) = (1, 10) at Re = 75 x 10"^, an arrangement which will be shown 
later to be representative of a VLSM. The normal wave speed, Uw = a;/(fc^ -I- n^)^/^, 
is increased from 0.2 to 0.6 from figure 4 to figure 6. Note that because (fc, n) is held 
constant between the figures the wavespccd in the streamwise direction, Ux, is a constant 
fraction of the normal wavespeed, and that the choice of fc 7^ 0, n 7^ leads to helical 
mode shapes. As expected, the energy for the lower wave speed is concentrated close to 
the wall, in the inner (wall) region, while for the faster wave the energy is centred in the 
core of the pipe. This trend is to be expected from Taylor's frozen turbulence hypothesis 
if propagating waves are a realistic feature of wall turbulence. 

For this wavenumber pair, the forcing elicits an energetic response that is concentrated 
in the streamwise component for each lu, with the relative magnitudes of the azimuthal 
and radial velocities being frequency-dependent. The modes for the lower two frequen- 
cies shown in figures 4 and 5 have footprints which reach down to the wall. Increasing w 
leads to a "lifting" of the streamwise and spanwise velocity components, manifested as 



A critical layer model for turbulent pipe flow 



13 




Figure 2. Radial variation of the absolute amplitude of each component of velocity for the 
response modes with {k,n,uj) = (2,2, 1) associated with the first three singular values, am, at 

Re = 75 X 10''. — : abs{u); : abs{v); ■ ■ ■ : abs{w). Top, middle and bottom panels correspond 

to m = 1, 2, 3, respectively. 




Figure 3. Variation of singular values am with mode number, m, for (fc,n) = (1,10) and 
Re = 75 X 10"^. For a wide range of frequencies, the first singular value is an order of magnitude 
larger than the second. 



a distinctive inclination to the wall of velocity isocontours. At the highest of the three 
frequencies, shown in figure 6, the mode detaches from the wall and the streamwise and 
spanwise velocities display something more like a two-level structure than inclined isocon- 
tours. For all three frequencies, the wall-normal component shows little phase variation 
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with wall-normal distance. The shapes imply a structure of long rolls with strcamwise 
vorticity and strong streaks in the streamwise velocity. Note that the velocity distribu- 
tions shown in figure 5 are in excellent agreement with the conditionally-averaged very 
large scale structures determined for channel flow by Chung & McKeon (2009) and the 
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streamwise coherence inferred in from two-point correlations with a reference point out- 
side the immediate near wall region in a high Reynolds number boundary layer by Guala 
et al. (20096), suggesting that modes of this kind may have particular importance for 
wall turbulence. Section 3.4 investigates modes with self-similar energy distributions in 
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Figure 6. Structure of the velocity field for the first mode with {k,n,ij) — (1,10,0.6) at 
Re = 75x 10^. Top panel s: u, V and w in the streamwise-wall-normal plane. Bottom: variation 
of modal kinetic energy. 



and y/R for various wavcnumbcr pairs. Wc call these "inner" and "outer" scaling 
modes respectively. 
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3.2. Effect of streamwise wavespeed on the radial distribution of perturbation energy 

In this section we explore the influence of frequency, or more accurately the streamwise 
wavespeed Ux/Ucl^ on the mode shapes. The radial distribution of energy in each ve- 
locity component for the mode with (fc,n) = (27ri?"'"/1000, 27ri?"'"/100) and streamwise 
wavespeed in the range 0.1 ^ Ux/Ucl ^ 1 at Re = 410 x 10'^ is shown in the composite 
contour plots of figure 7. This wavenumber combination, (Aj, A+) ~ (1000, 100) with an 
appropriate phase velocity ~ 10, appears to be representative of the signature of 
the near-wall cycle in many studies across a range of Reynolds numbers and flow types. 
In other words, it is a universal feature of wall turbulence responsible for the near- wall 
peak in the streamwise intensity and the wall shear stress. The changing amplitudes in 
each velocity component as the wavespeed increases reflect the distribution of energy 
between components and across the radius for different frequencies. The distributions 
are normalised to give an identical total kinetic energy for each mode. 

The figure shows that for low wavespeeds the contours of constant energy lie at ap- 
proximately the same wall-normal locations for all three velocity components and that 
the mode reaches down to the wall (i.e. is "attached"). The distribution of energy be- 
tween the velocity components is consistent with the well-known picture of streamwise 
rolls and streaks, in which the wall-normal and (double-peaked) azimuthal velocities 
have a similar order of magnitude and give rise to much larger fiuctuating streamwise 
velocities. This type of distribution is observed for this wavenumber combination at all 
Reynolds numbers considered, although high radial resolution is required to capture it 
accurately as the Reynolds number increases. This is directly related to the experimental 
and computational difficulties associated with an increasing range of scales. 

A transition occurs in the radial and component- wise energy distribution with increas- 
ing wavespeed. Instead of the near- wall roll/streak pattern described above, the peak 
energy moves away from the wall and the energy is concentrated in the core region. The 
wall-normal and streamwise energies arc of similar order of magnitude (larger than the 
azimuthal component of energy), such that the mode pattern is more indicative of the 
signature of concentrated spanwise vorticity. 

Aspects of the results for this wavenumber pair arc representative of the trends ob- 
served for other modes. The location of the peak energy associated with a particular 
mode moves slowly away from the wall as the frequency is increased. For a range of 
low uj {Ux), the energy is concentrated in the near-wall region. For higher uj, the energy 
is concentrated in the core of the pipe. The shape of the velocity profile near the wall 
constrains the range of possible mode shapes. Specifically, it constrains the radial extent 
of the mode, providing a minimum distance from the wall for the energetic peak in the 
streamwise stress. For example, in figure 7, this occurs at y+ ~ 10. The wall-normal 
distance to the peak is Reynolds number dependent. 

3.3. Comparison of streamwise wavespeeds and local mean velocity 

In the following discussion, we take the wall-normal location of the streamwise energy 
peak to be representative of a nominal mode centre. This allows comparison between the 
streamwise component of wavespeed and the local mean velocity. From figure 7, we expect 
that these will not always be equal. Figure 8 shows the variation of peak energy location 
with streamwise phase speed, for {k, n) = (1, 10) and Reynold numbers of 75 x 10'^ and 
410 X 10"^. Whilst there is a region in the core of the pipe where the wave has a streamwise 
wavespeed that is similar to the local mean velocity at the peak of the streamwise energy, 
close to the wall there is a significant deviation between the two velocities. This implies 
that our assumption that the mode is centred at the streamwise energy peak is incorrect 
near the wall, that these modes do not obey Taylor's hypothesis , or that the modes 
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Figure 8. Variation of location of peak streamwise perturbation amplitude with streamwise 

phase speed for modes with {k, n) = (1, 10). — : iie = 75 x 10^; Re = 410 x 10^. The dashed 

gray lines show the turbulent mean velocity profiles at each Reynolds number. Right panels: 
variation of resolvent norm with streamwise phase speed. 



will not be observed in the flow. In order to distinguish between the distinct behaviour 
near to and very far from the wall, we designate the former "wall modes" and the latter 
"critical modes" , for reasons that will become clearer in section 4. The right-hand panel 
of figure 8 shows that the first singular value for this wavenumber pair increases with 
Reynolds number, at least in the region far from the wall. 

Figure 9 shows the influence of azimuthal wavenumber on the streamwise wavcspeed 
corresponding to the peak in streamwise energy, for a mode representative of the near- wall 
cycle at Re = 75 x 10^ and another mode with {k,n) = (1, 10) at Re = 410 x 10^, re- 
spectively. For azimuthal wavenumbers that are not "too large" , both cases show critical 
modes in the core of the pipe over a range that increases with increasing n, with wall 
modes observed closer to the wall and a transition region in between. However, when the 
azimuthal wavenumber becomes sufficiently large, only wall modes are observed. Then, 
the peak energy remains at an approximately constant wall-normal location indepen- 
dent of wavespeed and the mode shape is self-similar with decreasing wavespeed. In this 
case, the corresponding singular value is small and almost independent of streamwise 
wavespeed. 

Wherever the streamwise wavespeed is equal to the prescribed local mean velocity, 
the behaviour approaches that of a "critical" layer. Very large amplitude, spatially and 
temporally periodic waves have been observed in various shear flows (Maslowe 1986) and 
explained in the context of critical layers. 

Examination of the resolvent (Equation 2.15) shows that the system response ap- 
proaches singular for a phase velocity close to the local mean {oj/k ^ U) and for high 
Reynolds number. That is, the flow approaches a state where it supports neutrally sta- 
ble modes, corresponding to an eigenvalue approaching the imaginary axis. Comparison 
between the right-hand panels of figure 8 shows that the resolvent norm is orders of 
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Figure 9. Variation of location of mode peak amplitude with spanwise wavenumber. Top: 

A+ w 1000 and = 100 (— ), 50 ( ), 30 (• • • ) and 10 (• - •) at i?e = 75 x 10^ Bottom: k = l 

and n = 10 (— ), 100 ( ), 200 (• ■ • ) and 400 ( ) at Re = 410 x 10^. The dashed gray lines 

show the turbulent mean velocity profiles at each Reynolds number. Right panels: variation of 
resolvent norm with streamwise phase speed. 

magnitude larger for such critical modes. This is phenomenon is explored further in 
Section 4.2 below. 

3.4. Inner and outer scaling modes and the effect of Reynolds number 

By carefully selecting the streamwise wavespeed, we observe scaling behaviour that is 
surprisingly consistent with that seen in experimental velocity spectra. The modes ap- 
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proach self-similarity when scaled on the classical inner or outer scales. Figure 10 shows 
the {u,v,w) mode shapes for a wavenumber-frequency combination that is representa- 
tive of the near-wall cycle, (k, n) w (27ri?+/1000, 27ri?+/100) with strcamwise wavcspeeds 
equal to = 10 and 20. The former velocity is known to correspond to the convection 
velocity of near-wall structure, while the latter velocity will only correspond to the in- 
ner scaling region for sufficiently high Reynolds number, when the latter velocity occurs 
below the outer edge of the log region in the mean velocity. It is understood that the 
velocity spectra at high wavenumbers near the wall collapse when scaled on inner units, 
while the lack of collapse of the intensities arises from the influence of larger scales (Met- 
zger & Klewicki 2001). Similarly, the intensities approach collapse in the core of the pipe 
as the Reynolds number is increased (McKeon & Morrison 2007). 

Figure 10 shows that these modes become self-similar when scaled in inner units across 
nearly two decades in Reynolds number and that they are "attached" in the sense that 
their footprints reach down to the wall. Remarkably, the peak energy in the strcamwise 
u component for the lower wavespeed of lOu^ is at the appropriate wall- normal location, 
y"*" « 15 — 20; obtained from experimental and computational observations of the near- 
wall cycle. The increase in the maximum power occurs because the modes are co-located 
in plus units, so their extent decreases in dimensional units as the Reynolds number 
increases, and the total perturbation energy in each mode is normalised to one. There is 
an upper limit to the range of Reynolds number that can be considered due to numerical 
resolution. This is discussed in Appendix B. However the Reynolds number trend appears 
to be quite clear. 

A similar result can be obtained for modes that occur in the core of the pipe, i.e. they 
are expected to scale on outer variables. Figure 11 shows the {u,v,w) mode shapes for 
(fc,n) = (1,10) and constant velocity defects, U^j^ — U^. The approach to similarity 
with increasing Reynolds number is slow, in agreement with the scaling of the integrated 
turbulence intensities (McKeon & Morrison 2007), with the wall- normal component col- 
lapsing earliest. 

We see that a straightforward analysis of the resolvent, a linear operator, produces 
mode shapes that agree with known scaling remarkably well. In one respect this simply 
reflects the self-similarity of the mean velocity profile near the wall, in the sense that 
the modes are a result of the local shape of the velocity profile, which does not vary in 
inner units once the Reynolds number is sufficiently high. However the recovery of this 
result underlines the utility of the model for describing turbulent fiow. A more formal 
justification for the scaling will be given in the discussion presented in Section 4. 

As evidenced by the failure of inner scaling to collapse experimental observations of 
the integrated strcamwise energy near the wall, the transition from inner to outer scaling 
is complex and the range of modes exhibiting inner and outer scaling at a given fiow 
condition is a function of the Reynolds number. We explore this regime further in the 
next section. 

3.5. The very large scale motions and transition from inner to outer scaling 

The existence of inner and outer scaling modes in our model is consistent with the classical 
turbulence scaling picture, in which viscous and outer scales are sufficient to describe the 
turbulent behaviour in the inner and outer regions of the flow. However, as described 
above, recent work suggests that there is some outer influence on the inner structure. For 
instance, Metzger & Klewicki (2001); Hutchins & Marusic (20076) show that the position 
of the near-wall peak in strcamwise energy weakly depends on Reynolds number. We 
consider here a possible origin for this effect, focusing on the transition between inner 
and outer scaling of the velocity modes. 
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Figure 10. Energy distribution over the pipe radius for {k,n) = (27r_R+/1000, 27rJ?+/100), 
or equivalently, {\t,^t) ~ (1000, 100). The equivalence is not exact due to the restriction of 
integer azimuthal wavenumber. The left panels show streamwise wavespeed — 10 and the 
right panels U+ = 20. Reynolds numbers: — 75 x 10^ 150 x 10^, • • • 410 x 10^, • - • 1 x 10**. 



A concatenation of experimental results in wall turbulence in general, and specifically 
in pipe flow, have shown that a streamwise wavenumber A: = 1 is representative of the 
VLSM phenomenon (Kim & Adrian 1999; Morrison et al. 2004; McKeon & Morrison 
2007; Hutchins & Marusic 2007a; Monty et al. 2007) although the recent study of Monty 
et al. (2009) suggests that the exact details vary from flow to flow. The appropriate 
spanwise wavenumber is less clear; the spanwisc measurements in pipe flow of Monty 
et al. (2007) suggest « 0.5R (n ^ 12), while the POD analysis and correlations of 
Bailey et al. (2008) indicate that lower values, n ^ 3, arc more realistic for the VLSMs. 
We investigate n = 10 as a compromise that follows the expected critical layer aspect 
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Figure 11. Energy distribution over the pipe radius for (fc, n) = (1, 10). The left panel shows 
the distribution for streamwise wavespeeds such that the velocity deficit is i7^^ — = 4 and 

the right panel for a velocity deficit of (7^^ — = 2. Reynolds numbers: — 75 x 10'^, 

150 X 10^ • • • 410 X 10^, • - • 1 X 10'^. 



ratio discussed in section 4 and is in reasonable agreement with the work of both Monty 
and Bailey. 

The variation of mode shape for {k, n) = (1, 10) with increasing phase speed is shown 
in figure 12 for Re = 410 x 10"^. We have defined a mode as critical when the streamwise 
wavespeed reaches the local mean velocity at the peak modal energy. The variation of 
this position with Reynolds number for the wavenumber pair related to the VLSM is 
illustrated in figure 8. In this case the mode first becomes critical when Ux is a con- 
stant fraction of the centreline velocity, namely Ux/Ucl ~ 2/3, independent of Reynolds 
number. This relationship suggests that for the critical modes the appropriate scaling ve- 
locity is the centreline velocity, rather than the friction velocity. This was also proposed 
by, among others, Jimenez et al. (2004), who observed that the so-called "global" modes 
in channel flow simulations appeared to convect with a velocity equal to 0.5C/cl- 

The foregoing section has demonstrated that the simple model for pipe flow is capa- 
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ble of reproducing several features of wall turbulence, from a self-similar distribution 
of energy for small scales near the wall to a large scale modal shape that is reminis- 
cent of recently-observed conditional averaged structure in channels (Chung & McKeon 
2009), two-point correlations in medium and high Reynolds number boundary layers 
(Guala et al. 20096), and energetic POD modes in pipe flow (Duggleby et al. 2007). In 
a broad sense, the singular value decomposition selects the mode shape at each (fc, n, uj) 
combination that is "most likely" to be observed based on receptivity to forcing at that 
combination, or "optimal" in a sense analogous to that of earlier initial value studies. The 
amplitude of the response is set by the amplitude of the forcing, which is determined by 
other (fc, n, uj) combinations such that the ensemble of mode shapes, and specifically the 
Reynolds stress distribution, is consistent with the base mean flow. Note that a canonical 
spectral analysis records the integral energy over all modes with the same wavenumbers, 
so there is not a one-to-one relationship between the information given by the different 
bases corresponding to the power spectrum and the modal decomposition, except in an 
integral sense. 



4. A wall and critical layer framework 

In this section, we show that the simple model presented in Section 2 predicts the 
location of the peak in streamwise turbulent energy associated with a wavenumber pair 
representative of a VLSM-typc motion. We then offer a critical-layer interpretation of 
this type of mode. 

4.1. Scaling of the VLSM energy peak 

The peak energy in the streamwise velocity component for the first critical mode with 
(fc,n) = (1, 10) occurs for y/R < 0.1 for all the Reynolds numbers considered here. Thus 
it is clear that the first critical mode occurs within the log region, at least for sufficiently 
high Reynolds numbers. f Based on this velocity scaling, the location of the energy peak 
associated with this mode can be predicted using similarity of the mean velocity in both 
the inner and core regions, in the form of a logarithmic profile in the overlap region and 
Reynolds similarity of the outer flow. The inner-scaled deviation of the centreline velocity 
from a log law is a constant, C. The relationship is 

f/+(2/+) = -lny+ + S, (4.1) 

K 

U+^^U+{y+ = R+) = -\nR+ + B + C. (4.2) 

The velocity reaches a value of two-thirds of the centreline velocity at an inner-scaled 
wall-normal distance yy^ that can be predicted as follows: 

U+ = - hiy+ +B = l(-\nR+ + B + c] , (4.3) 

or, rearranging, 

y+3 = ai?+2/3 (4.4) 

where a = exp[f (2C - B)]. For a pipe k = 0.421, B 5.60 (McKcon et al. 2004) and 
C w 2 (McKeon et al. 2005), so a = 0.8. 

t A logarithmic profile will give a reasonable approximation for the local velocity at the 
lower Reynolds numbers, even if the profile does not correspond to an overlap region that is 
independent of Reynolds number. 
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Figure 13. Comparison of experimental pipe flow results on the y -location of the spectral 
peak in streamwise energy associated with the VLSM modes and the prediction of equation 4.4. 
• : Superpipe data (Morrison et al. , 2004); dashed Une: i/+ = 0.8-R"'"^''^. Note that the high- 
est Reynolds number data point may have been affected by hot-wire spatial resolution (the 
non-dimensional wire length at this Reynolds number is l'^ — 385). 

The agreement is excellent with both the (fe, n) = (1, 10) mode peak and the previously- 
pubhshed experimental variation of the VLSM streamwise energy peak shown in figure 13. 
For the experimental data, the location of the VLSM peak at each Reynolds number was 
determined by considering the low wavenumber spectral peak in the hot-wire data set of 
Morrison et al. (2004), as reported in McKeon (2008). It should be noted that identifying 
the location of the peak energy relies on the local smoothness of the spectra. Thus the 
error bars on the peak position shown in the figure arc conservative and correspond to the 
next closest wall-normal locations at which spectra were obtained. While probe resolution 
effects may be a concern at the highest Reynolds number, the data at i?"*" = 19 x 10^ 
clearly represent an outlier. For the remaining data the agreement is excellent. Thus we 
associate the first critical mode at this wavenumber pair, namely the critical mode with 
the lowest phase velocity, {k,n,u!) = (1, 10, 2/3J7cl), with the experimentally-observed 
VLSM, and infer that this mode gives a dominant energetic contribution to the turbulent 
fluctuations. Equation 4.4 represents a non-observational attempt to predict the location 
of the peak in the streamwise energy associated with the VLSMs. 

The scaling with the two-thirds power of Reynolds number is also reminiscent of critical 
layer arguments for neutrally-stable disturbances in linearly unstable laminar flow, thus 
we expand the critical layer framework to include forced, propagating modes in turbulent 
pipe flow. 

4.2. A wall and critical layer framework 

A consequence of the decomposition performed in Section 2 is that equation 2.11 can be 
rewritten in terms of the vertical velocity and vorticity, yielding an operator equivalent 
to the more familiar Orr-Sommerfeld-Squirc operator when the flow is linearised around 
the turbulent base flow, identified earlier as the {k,n,u]) = (0,0,0) mode. As such, we 
extend the tools of linear stability analysis to explore the scaling of the layers around the 
critical points, with reference to known results concerning the OSS equations. 

To simplify the exposition, we consider the case for plane flow in Cartesian coordinates. 
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Linearised plane flow has unstable (right-half plane) eigenvalues for sufficiently high 
Reynolds number and for the unforced case the "neutral curve" bounds the region of 
stability in the Re — k plane. Pipe flow is always linearly stable, so has no neutral curve. 
However, comparable physical processes occur and are manifested as regions of high 
resolvent norm without actually reaching linear instability. This said, we will work with 
the standard formulation of the Orr-Sommerfeld equation (0-S) for a plane flow to aid 
the discussion, 

{U - a;/fc)(p2 _ K^)y _ u"v - -^{V^ - K^v = 0. (4.5) 

ikRe 

Here U" is the second derivative of the base profile. 

Scaling analysis of this equation reveals two wall-normal regions in which the action of 
viscosity is required. The singularity in the inviscid 0-S (Rayleigh) formulation occurring 
if Lu/k ~ U may be resolved by restoring either viscous effects or nonlinearities in a 
region local to the critical point. We consider the case of viscous critical layers here, 
with justification to follow. Since the solutions of the Rayleigh equation do not obey the 
viscous boundary conditions, a second region in which the effects of viscosity are restored 
is required close to the wall. The scaling of the width of these regions can be determined 
by consideration of the appropriate terms in Equation 4.5 (for plane flow see Schmid 
(2001); Drazin & Reid (2004); Maslowe (1981)). 

Close to the wall, the boundary layer approximation in a region around yi is given by 

-^V^i = -uj/kV^i! = -U^V^v 
ikRe 

with solution 

i{y) = Cexp (-(y - yi){ikReU,)^/^e-'^/^') (4.6) 

such that the viscous layer around yi has thickness of order kRe^/'^. 

At the critical layer centred on yc the approximation to the 0-S equation is 

which reduces to an Airy equation for V" 

under the substitutions e = {ikReU[)~^/'^ , £, = [v ^ Vcilf- and v(y) ~ y{C)- Therefore 
the critical layer thickness is of order yjR ^ kRe~^l^ and the mode phase velocity, 
Vx = w/fc, sets the wall- normal position of yc- 

Thus we expect to observe viscous and critical modes with (kRe~^l^) and {kRe~^l'^) 
scaling, respectively, in the wall-normal direction. Note that in inner scaling, this corre- 
sponds to 2/+ ~ (i?+/fc)-^/^ and 2/+ ~ (i?+^/fc)^/^, respectively. 

The two options for the relative locations of the two layers are shown schematically in 
figure 14: if the layers are distinct (shown in the left panel) then the critical layer exists 
at a wall-normal distance 0(kRe)~^l^ and if the critical layer reaches down to the wall, 
then it will be centred at a wall-normal location 0(kRe)~^l^ (shown in the right panel). 
In the 0-S analysis, these two cases correspond to the upper and lower branches of the 
neutral stability curve at the Reynolds number under consideration. 

A similar scaling analysis can be performed in turbulent pipe flow using the turbulent 
base profile and with an appropriate Reynolds number for turbulent flow. Re = i?"*". 
There are no neutral curves because the eigenvalues remain in the left half plane for all 
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Figure 14. Schematic of the scahng of upper (left) and lower (right) branch critical and wall 

layers in the Orr-Sommerfeld equation. 

Reynolds number and all wavenumber-frequency triplets. However similar effects to the 
true 0-S system are manifested as high system response to forcing, giving a high resolvent 
norm with the same scaling. Hence we refer to upper and lower branch-type modes in 
what follows and consider the implications of wall and critical layers for the scaling of 
turbulent pipe flow and in particular the implied mechanism for the manifestation of 
viscous effects outside of the immediate near-wall region. 

We conclude the this part of the discussion by observing that in laminar pipe flow 
the Reynolds-number independence of the form of the velocity profile means that there 
is a linear relationship between the point at which the velocity reaches two-thirds of 
the centreline value and the Reynolds number, ~ . It is known that the linear 
Navier-Stokes operator is particularly sensitive to perturbation at this point, where the 
wall and centre mode branches of the spectrum merge (Rcddy et al. 1993). 

4.3. Non-normality versus criticality 

The foregoing analysis permits a comment on the importance of non-normality to the 
sensitivity to forcing. As indicated in Section 2, the resolvent in Equation 2.15 can become 
large in two distinct ways, via the non-normality arising from the gradient of the mean 
velocity, or when the diagonal terms tend to zero at a critical layer. Regarding the 
non-normal mechanism, the local shear couples the streamwise and wall normal velocity 
components, leading to a "lifting" mechanism for streamwise vorticity. On the other 
hand, the critical layer arguments arise from arguments that are normal in character, 
since this mechanism is present even when the mean shear is small. This provides a simple 
explanation of why critical modes are observed far from the wall in our model and become 
increasingly dominant at higher Reynolds number. The decomposition performed in this 
work therefore permits a simple identification of the two effects. 

4.4. VLSMs as a consequence of critical layer scaling 

We predicted in Section 4.1 that the wall- normal location of the peak energy of the 
mode associated with the very large scale motions scales with R^^/^. The component- 
wise velocity distributions for the VLSM mode have been shown in figure 5. Many of 
the characteristics of this mode are in good agreement with those of the lower branch- 
type critical layer under the analysis described above. In turbulence terminology, the 
associated disturbance would be "attached" to the wall in the sense that it has finite 
amplitude in the near- wall region. This property is consistent with the observations of the 
"footprint" of the VLSM modes reaching down to the wall (Hutchins & Marusic (2007a), 
Guala et al. (20096) and others) and influencing the instantaneous wall shear stress 
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(Marusic & Heuer 2007) . This is a consequence of this mode retaining characteristics of 
the wall layer while having the lowest phase velocity which can be considered critical. 
This phase velocity is a constant fraction of the centreline velocity, in agreement with the 
conclusions of del Alamo et al. that the centreline velocity is the appropriate velocity 
scale for the very long, "global" modes. 

The mode shapes associated with the VLSMs in figure 5 are in good agreement with 
the structure that can be inferred to give rise to the two-point correlations reported in the 
near-neutral surface layer by Guala et al. (2009 &) and the conditional averages of Chung 
& McKeon (2009) in channel fiow. This suggests that the VLSM mode predicted here 
indeed becomes a dominant, viscous structure in the near-wall region as the Reynolds 
number increases. 

The VLSM mode predicted here also has a concentration of fluctuation energy in the 
wall-parallel components, in agreement with the experimental evidence that the very 
long scales arc visible in the u and w spectra (although the literature on the latter is 
relatively scarce due to the experimental difficulty of obtaining well spatially resolved 
information on the spanwisc component), but not in the wall- normal component due 
to the effect of "blocking" due to the wall on eddy distribution. In addition, the spatial 
variations of the streamwise and wall-normal velocities implies a distribution of Reynolds 
shear stress that demonstrates the tt phase reversal associated with viscous critical layers 
in the vicinity of the peak in the streamwise energy. This justifies our earlier assumption 
that viscosity is more important than nonlinearity in the vicinity of the critical layer. 
The phase relationships between the axial and wall-normal velocity components implies 
that these scales are "active" (in the sense of Townsend (1976)) in that they bear non- 
negligible Reynolds shear stress, as proposed by Guala et al. (2006), but in contrast to 
the spirit of Townsend's original ideas. 

We note also that counter-rotating vortex structure implied by the VLSM mode shape 
(figure 5) is a well-known phenomenon associated with lower branch critical layer solu- 
tions, at least in the laminar case (Viswanath 2009; Wang et al. 2007), and are observed 
at similar wavenumber pairs, n ^ lOfc, using the current model for laminar fiow (Sharma 
& McKeon 2009). The conditionally-averaged cross-stream streamlines and swirl distri- 
butions of Hutchins & Marusic (2007&) and Chung & McKeon (2009) confirm that this is 
the expected mode shape for the VLSMs in channel flows. This structure is also consistent 
with the work of del Alamo & Jimenez, who showed that the two types of disturbance 
that experience maximal transient growth in the initial value problem in turbulent chan- 
nel flow, one inner scaling and one outer scaling mode, both resemble streamwise rolls 
and streaks. 

We conclude our comments on the relationship between VLSMs and critical layers 
by proposing a partial explanation of the apparent meandering, very long streamwise 
coherence of VLSMs that is indicated by experimental measurements of streamwise ve- 
locity in the spanwise plane for various flows. To this end, figure 15 shows isosurfaces of 
streamwise velocity for a sum of left- and right-going propagating VLSM-type modes. It 
seems plausible that the superposition of other modes on such a pair would reproduce 
something like the apparent observed meandering coherence of the order of 25(5 (Monty 
et al. 2007). 

4.5. The importance of wall layers 

Wall layer scaling has been implicit in the work of several researchers, including the 
experimental and theoretical discourses by Sreenivasan & Sahay (1997) and Sahay & 
Sreenivasan (1999) and the form of the mean velocity in the vicinity of this peak proposed 
by Sreenivasan & Bcrshadskii (2006). The proposed mean momentum balance structure 
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Figure 15. Isocontours of streamwise velocity for the sum of left- and right-going VLSM 
propagating modes, {k,n,uj) = (1, ±10, 1 * 2/3Ucl)- The Reynolds number is Re = 75 x 10^. 



of Klewicki et al. (2007) describes the scaling in inner units with of the Reynolds 

stress peak location. An analogy between the wall layer scaling phenomenon and linear 
stability concepts was drawn by Sreenivasan (1988) who explained the dependence 
in terms of the roll-up of a vortex sheet used to model the vorticity in the boundary 
layer. Sirovich et al. (1990) observed plane waves in channel flow travelling at the speed 
of corresponding to the mean velocity at the location of the Reynolds stress peak. 

The foregoing discussion has given a rigorous derivation of the origin of the wall layer 
scaling. Such arguments can be extended to discuss local scaling of the turbulent spec- 
trum. While it is understood that the turbulent fluctuations in the near wall region resist 
inner scaling, the cause of this discrepancy has been understood for some time to lie in 
the increase in energy at larger streamwise wavenumbers as the Reynolds number is in- 
creased, with the implication that the spectral shape at large k (short wavelengths) has 
the potential to be self-similar. In fact, the discrepancy is a consequence of wall layer 
scaling, which we now demonstrate. 

Applying wall layer scaling for a mode with streamwise wavenumber k gives 

For a constant A J , this can be rewritten in the following way, demonstrating that the peak 
energy location scaled in inner units will be independent of Reynolds number, provided 
that the mode phase velocity also remains constant: 

1/2 

(4.9) 



Vpk 
R 



1/2 



2tt 



Based on the results shown earlier for A+ sa 1000, where y^f, w 13, the constant of 
proportionality is 0(1). So for the wavenumber range corresponding solely to wall layer 
arguments, the spectral shape should be self-similar in inner scaling. This was confirmed 
in figure 10. Predicting exact self-similarity would require further arguments relating to 
the Reynolds number dependence of the product of the forcing amplitude and the first 
singular value, which is beyond the scope of this paper. 
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4.6. Amplitude modulation and the phase relationship between the large and small scales 

A region in the flow where an apparent amplitude modulation of the small scales occurs 
has been identified and explored by Blackweldcr & Kovasznay (1972); Bandyopadhyay 
& Hussain (1984); Hutchins & Marusic (20076); Mathis et al. (2009) and Guala et al. 
(20096). This region can be identified with the concatenation of viscous layers at different 
wavenumber/frequency combinations and is required to meet the wall boundary condi- 
tions. Thus it must respond to the influence of finite amplitude critical layer disturbances. 
It may be said, then, that it truly represents an amplitude modulation of turbulence ac- 
tivity at the range of scales present in the near wall flow, the range of which is known 
to correspond to wavelengths smaller than the VLSMs. By the arguments made above, 
this viscous layer, in which the modulation is largest, should scale with This is 

borne out in the near-neutral atmospheric surface layer with (5+ ^ 0(10^) in the mea- 
surements of Guala et al. who show that the difference in energetic content of the small 
scales between positive and negative excursions of the fluctuation velocity in the VLSM 
bandpass-filtered signal is concentrated in the region y"*" < O(IO^). In addition, the cor- 
relation coefficient between the large scale component and the filtered envelope of the 
small scale component of Mathis et al. (2009) shows a distinct discontinuity in gradient 
that tracks a position y+ ~ 5+^/^. We hypothesise that this corresponds to the outer 
edge of the inner viscous layer associated with the most energetic critical layer mode, 
which we identify as a VLSM. 

4.7. Implications for the scaling of the fluctuations in turbulent pipe flow 

Previous phenomenological observations and the predictions of the critical layer scaling 
can be concatenated into a skeletal description of turbulent pipe flow. We showed above 
that wall layer scaling implies a similarity over the part of the turbulent spectrum gov- 
erned by wall modes. Additional arguments for the outer region of the flow (in which only 
critical modes contribute turbulent energy) seem to require the exact form of the mean 
velocity profile. However the similarity of the mode shapes for a constant velocity defect 
at high Reynolds number in figure 11 suggests that outer scaling of the spectrum can 
be recovered. Therefore while outer fluctuations should be close to self-similar in outer 
units, the inner fluctuations will always include contributions from both wall and critical 
layer modes. This, perhaps, is one cause of the apparent success of the "mixed scaling" 
proposed by DeGraaff & Eaton (2000), which leads to collapse of m^/ [urU oo) over a lim- 
ited Reynolds number range. In light of the distinct footprints of the different classes of 
modes, the controversial question of the importance of "top-down" versus "bottom-up" 
effects may be more appropriately framed in terms of the relative importance of wall and 
critical layers, with the former tied to the latter through the wall boundary conditions, 
but with all fluctuations contributing to the shape of the mean profile. 

The distributions of Reynolds shear stress associated with wall and critical modes have 
distinct shapes and, specifically, are local in That is, the Reynolds shear stress is also 
located in the perturbation layer, close to the peak energy. In addition, this localisation 
goes towards explaining some of the observations of Guala et al. (2006) of the VLSMs, 
that these low wavenumbers (associated with critical modes in our formulation) are active 
in the sense that they contribute significantly to the Reynolds shear stress, and also that 
the sign of the wall-normal gradient in shear stress at this wavenumber, A[—uv{k))/Ay, is 
opposite to that of the smaller scales closer to the wall. The reversal in sign of the mean 
Reynolds stress associated with a single mode lends support to the notion of a spatial 
transfer of turbulent energy away from the wall to the outer region, by simple arguments 
concerning the sign of the local turbulent energy production. Of course, in the overlap 
region the integrated shear stress over all modes is such that local equilibrium concepts 
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are at least approximately met. There is a coupling between the sum of Reynolds shear 
stress contributions from all turbulent modes and the mean velocity profile that yields 
those modes. 

For the VLSM mode, at least, a logarithmic variation in the mean velocity implies 
that this critical layer is always associated with a constant ratio of local to centreline 
velocity, {U/Ucl)~^, with the consequence that knowledge of the location of the centre of 
the VLSM mode, which can be experimentally determined from the peak in the stream- 
wise turbulent energy, yields Equation 4.4 relating the log law constants k, B and C. 
Alternatively, if the log law constants arc known, then the peak offers an alternative way 
to determine the friction velocity, Ur- While the footprint of the VLSM in the streamwise 
spectrum becomes more dominant as the Reynolds number increases, accurate identifi- 
cation of the location of the maximum becomes more difficult because the distribution 
around the peak also broadens in y"*", so these relationships may not prove to be useful 
in a predictive sense. 

Our key assumption is that the modes indicated by the singular value decomposition 
will be observed in turbulent pipe flow if the wall-normal distribution of the forcing at 
that (fc,n,cj) combination is non-zero, contains an exactly appropriate component and 
the singular value is large. Clearly there are many other modes with the same streamwise 
wavenumbers or frequencies, as observed in the spatio-temporal spectra of, for example, 
Morrison & Kronauer (1969), but different convective velocities which blur the harmonic 
decomposition. However, Hutchins & Marusic (2007a) have used a synthetic signal to 
show that an isolated portion of a signal with clear streamwise and spanwise spatial 
content leads to a broad peak in the spatial power spectrum centred on the true signal 
wavenumbers. We propose that the other modes combine to mask the signature of the 
VLSMs when a Fourier decomposition is performed, such that the mode corresponding to 
the dominant VLSM critical mode is not clearly deflned at fc = 1 in the power spectrum. 

We do not address the spatial phase relationship between modes, however we suggest 
that the base of knowledge concerning inclined ramp-like structures and the wall-normal 
arrangement of statistically-representative hairpin vortices, together with recent work on 
the interaction between the large and small scales (Mathis et al. 2009; Chung & McKcon 
2009; Guala et al. 20096) will inform future work in reassembling the modes predicted 
here to represent the essence of wall turbulence. 

4.8. Extension to other flows 

We have limited the arguments made above to turbulence in pipe flow. However an 
equivalent formulation to that detailed in Section 3.5 can be made for other flows by 
considering the best values of the Karman and additive constants, and the wake deviation, 
for other canonical flows and retaining the assumption that the critical layer occurs when 
U{y) = 2/3UcL (i-C. that the critical layer scales with y+ ^ For channel flow 

with K ^ 0.38, B = 4.17 and C = 0.175 (Nagib & Chauhan 2008), the estimation of 
the location of the peak associated with a lower branch-type critical VLSM is given by 
Equation 4.4, giving 

y+/3 « 0.62i?+2/3. (4.10) 

The difference in the constant arises because of the form of the mean profile close to 
the centreline (the wake component), rather than the subtlety of the log law scaling, the 
potential (lack of) universality of which remains under debate (Nagib & Chauhan 2008). 
The selection of A: = 1 to represent the VLSM structure in channel flow is supported by 
the work of Monty et al. (2009), but that there is also signiflcant energy observed at a 
streamwise wavelength of three channel half-heights, particularly in DNS studies. 
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Figure 16. Comparison of the {/"'"-location of VLSM-related phenomena in zero pressure gradi- 
ent turbulent boundary layers over a range of Reynolds number. •: location of zero amplitude 
modulation in laboratory boundary layers, estimates of the zero modulation location □ (both 
Mathis et al. , 2009) and the VLSM energetic peak ■ (Guala et al. , in preparation) with various 

predicted forms. — \/T5(5"'"^''^ (Mathis et al. , also wall layer prediction), 0.67(5"'"^''^ (lower 

branch-type critical layer prediction), ■ • ■ 0.765^*''^ (upper branch-type critical layer prediction). 



A similar value for the constant, a = 0.67 , is obtained by consideration of the log law 
constants in zero pressure gradient turbulent boundary layers. The analysis for this flow 
requires accounting for, or formally proving negligible, the influence of spatial inhomo- 
geneity in the streamwise direction, i.e. there is an implicit parallel flow approximation 
here. However there is utility in comparing the prediction of Equation 4.10 with the 
experimental data. The results are plotted in figure 16 which shows the experimentally- 
determined locations where the degree of amplitude modulation is zero in the data of 
Mathis et al. (2009) (which they equate with the VLSM energy peak) and several pre- 
dicted Reynolds number variations. Mathis et al. proposed S~^^^^ dependence of the peak 
location, after assuming that the VLSMs inhabit the exact middle of the overlap layer. 
We have shown earlier that it is also a consequence of wall layer scaling. The constants for 
the S+^/^ scaling were determined by consideration of the point where the local velocity is 
four-fifths of the freestream value, in a similar analysis to Equation 4.4. While the 5+^^^ 
and (5+^/'^ scalings could be said to describe the lower Reynolds number data equally well, 
only the S'^^^^ curve approaches one of the high Reynolds number data points obtained 
in the near-neutral atmospheric surface layer by Mathis et al. and Metzger & McKeon. 
The former study determined the higher point from a study with limited resolution close 
to the wall, while the latter study was performed with a maximum measurement height 
of the order of 10% of the equivalent boundary layer thickness (and can be seen to reveal 
two spectral peaks). The critical layer framework provides an interesting interpretation 
of these two identified spectral energy maxima. 

The data support a scaling with (J"*"^/^, particularly if the lower estimate of the VLSM 
peak energy in the surface layer is taken into account, suggesting that in the boundary 
layer there is a strong energetic contribution from the viscous wall mode (or, of course, the 
critical layer scaling changes due to the spatial inhomogeneity over the spatial VLSM 
period). The S~^^^^ scaling is in agreement with the observation by Jimenez that the 
global modes travel with a convective velocity Us ~ 0.5Ucl if the mean velocity is given 
by a log law governed by Equation 4.4. 

However the error bars around the peak in the streamwise turbulent energy identified 
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by Marusic et al. encompass the point y'^ = 0.76(5^'*/^, raising the intriguing possibility 
that we are observing the upper branch-type configuration of figure 14, witli two local 
maxima in the streamwise spectrum at locations proportional to S^^^"^ and 5+^/^. 
This may be a difficult hypothesis to prove, given the difficulty of obtaining fully-resolved 
spectral information for Reynolds numbers 20 x 10^ < (5+ < 7 x 10^ or distinguishing 
between two spectral peaks that are close in physical and spectral space. However we 
note that the streamwise velocity probability density functions reported by, amongst 
others, Adrian et al. (2000) and Morris et al. (2007) suggest two maxima that can be 
shown to correspond to approximately 0.5 and 0.8 times the free-stream velocity, which 
may confirm our hypothesis concerning the special significance of these velocities in the 
boundary layer. 

The existence of upper and lower branch-type modes in the different canonical config- 
urations of wall turbulence may contribute to the explanation for the emerging, subtle 
differences between the flows. In particular we note that the dominance of a different 
branch could affect the Reynolds stress profile with potential implications for the uni- 
versality of even the mean velocity profile. 

These results have assumed that the VLSM-type modes occur universally at fc = 1, 
while there is experimental evidence to suggest that the exact streamwise length varies 
from flow to flow, with the largest and smallest wavelengths observed in pipe and bound- 
ary layer flows, respectively. However the framework presented above can easily be mod- 
ified to account for this, should it prove to be required. 

5. Broader implications 

The critical layer arguments made here have several broader implications for canonical 
and more general flows. Next we identify some that are most pertinent to our current 
areas of study. 

5.1. Importance of the linear operator 

We have shown that consideration of the linear operator leads to a prediction for the 
dominant energetic mode shapes. Note, however that a linearisation per se is not strictly 
required. This result is in agreement with the understanding that the linear operator 
alone is responsible for the extraction of energy from the mean flow. This suggests that 
the modes predicted by this analysis correspond more to the production and spectral 
content of energy at a particular wavenumber pair, with the full spectrum determining 
the rate of transfer between wavenumbers via nonlinear interactions. 

5.2. On Taylor's hypothesis 

The modal decomposition proposed in this work is consistent with the limited obser- 
vations of the spatio-temporal spectra in the literature. The conversion between space 
and time is usually performed using Taylor's hypothesis of frozen turbulence, which is 
traditionally stated in terms of the slow evolution of the small-scale turbulence relative 
to the convective timescale and is believed to hold using the local mean velocity as the 
convection velocity for all but (possibly) the longest lengthscales (Dennis & Nickels 2008) 
and very close to the wall (Morrison & Kronauer 1969). The foregoing analysis suggests 
an alternative interpretation: the critical layer analysis implies that the dominant modes 
(and, by implication, structures) travel with a phase velocity that corresponds to the 
local mean at the peak streamwise mode energy. The reduction in amplitude away from 
the peak means that in the core part of the flow, the (fc — uj) joint spectra should be 
narrow with a peak ridge lying on a; = U{y)k. Nearer the wall, the influence of the wall 
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modes and the critical modes that extend to the waU but have energetic peaks further 
out will be to broaden the distribution of energy in (fc — w) space, remove the symmetry 
when integrating in k or u and ultimately lead to an inferred convcctive velocity that is 
different from the local mean velocity. 

This suggests that issues related to the validity of Taylor's hypothesis would be better 
expressed in terms of the symmetry and width of the k — cj spectrum. The work of Morri- 
son et al. (Morrison & Kronauer 1969; Morrison et al. 1971), among others, illuminated 
the deviation of the inferred convective velocity from the local mean, the lack of sym- 
metry near the wall in pipe flow and the spectral broadening that our model suggests. 
Recent Large Eddy Simulations in turbulent channel flow (Chung & McKeon 2009) reveal 
similar phenomena; this is an area of current experimental investigation. 

Consideration of the variation of the spanwise wavenumber of energetic modes, the 
n — uj spectrum, suggests an additional source for the spatial resolution problems that 
hamper the majority of laboratory high Reynolds number wall turbulence experiments, 
namely the spanwise extent of the sensing element in viscous units. 

5.3. The minimal unit for turbulence 

The "autonomous" simulations of Jimenez et al. (2004) demonstrated that a DNS of tur- 
bulent channel flow with a filter at a certain wall-normal height was capable of capturing 
the majority of the near-wall statistical features. This informed the "top-down" versus 
"bottom-up" debate and appeared to indicate that the external flow did not exert any 
control on the near-wall cycle. In light of the critical layer analysis, we observe that the 
wall-normal cut-off in these simulations, ^ 70, is above the peak of the VLSM critical 
mode for the nominal Reynolds number of Re-r ~ 550, for which ^ 50, so that even 
the lowest filter does not exclude the dominant, VLSM critical mode. To this extent, 
the results are consistent with our hypothesis. The implication for future simulations, 
then, is that the criterion on the minimum size of an autonomous simulation must be 
large enough to capture the azimuthal/spanwise and radial/ wall- normal extent of the 
critical mode associated with the VLSM, or the furthest reaching mode that leads to 
considerable modulation in a viscous layer near the wall. 

5.4. Experimental development length 

A simple estimate for the development length required in canonical wall turbulence ex- 
periments can be made by assuming the importance, or at least the large singular values, 
associated with modes with a streamwise-spanwise aspect ratio of 10:1, as in the VLSMs 
explored here. The peculiar constraint on azimuthal symmetry in pipe flow requires that 
n is integer, such that fc ^ 0.1 for modes with the required aspect ratio. This is in good 
agreement with the observed low wavenumber end of the premultiplied spectra in the 
overlap region in pipe flow (McKeon & Morrison 2007), at a wavelength A^^ « QOR. 
Combining this figure with the estimates detailed in Zagarola & Smits (1998) to account 
for the transition length at {L/R < 14) and the length for the shear layers to meet 
{L/R ^ 60) at Rejj ~ 10^ gives a reasonable estimate for the minimum required de- 
velopment length, Ld/R > 134_R. This figure is significantly exceeded in the Princeton 
Superpipe, where all measurements were taken for Ld/R > 164i?. The estimate is in 
good agreement with the length Ld ^ 160i? for fully-developed statistics in pipe flow 
(Monty et al. 2001) or i/i? w 134 in channel flow (Dean & Bradshaw 1976), both at simi- 
lar Reynolds numbers. This estimate also has implications for simulation, where box sizes 
can be set without consideration of the same entry effects experienced in experiments. 
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5.5. Turbulence over rough walls 

The wall and critical layers indicated by our model suggest a physical explanation for 
several keystones of our understanding of rough wall flows. Firstly, for flows that obey 
Townsend's hypothesis in the mean velocity, where the effect of the roughness is felt 
only in setting the friction velocity such that the self-similar form of the smooth wall 
is retained in the outer region, it is logical that the outer fluctuations would also obey 
similar scaling to the smooth wall case. This is because the outer region is dominated by 
modes whose shape is dictated by the local (i.e. outer) mean velocity profile. 

Secondly, while it is known that the near-wall streamwise intensity peak is suppressed 
relative to smooth wall values, the modulation effect associated with the VLSMs is still 
strong over walls with v/ur <?C fc < (5, for example in the near-neutral surface layer under 
conditions with k+ < 50 (Guala et al. 20096; Mathis et al. 2009). In this flow, the VLSM 
critical layer is centred outside of the roughness sublayer of order 3fc+ . The difference in 
the wall boundary condition would suggest that the structure of the viscous layer required 
to meet the wall boundary conditions would be quite different from that observed in the 
smooth wall case, in that it would be significantly less coherent and a strong function of 
the roughness geometry. However an analysis similar to that given in section 4.2 could be 
performed on the new mean flow, still yielding a wall layer that must meet the boundary 
condition in response to the velocity field imposed by critical modes. This is equivalent 
to stating that modification to the linear operator due to the roughness is confined to 
the region close to the wall, with the consequence that only the wall layer is significantly 
affected, for small enough roughness. Thus our critical layer model has the potential to 
make further progress in understanding the influence of surface roughness. 

5.6. Additional observations 

A full investigation of spectral scaling predictions using the critical layer model is re- 
served for future work. It also remains to cement the link between our model, statistical 
experimental data and the emerging picture of the structure of near wall turbulence. 
The picture has elements ranging from the VLSM structures discussed above through 
to observations of hairpin-like vortices, packets and uniform momentum zones (Adrian 
2007) and the near-wall cycle (previously thought to be autonomous) that is consistent 
with the physical arguments of Schoppa & Hussain (2002). Wc also note that the critical 
layer arguments hold significance for approaches to control of wall turbulence, specifically 
suggesting the possible utility of using distributed (harmonic) forcing at the wall. 

6. Conclusions 

We have described here a resolvent norm model capable of describing the most ampli- 
fied velocity modes in turbulent pipe flow, and an example of its application to illumi- 
nating scaling relationships in wall turbulence. 

To derive the model, we treat the nonlinearity in the perturbation equation (involving 
the Reynolds stress) as an unknown exogenous forcing, yielding a linear relationship 
given by the resolvent, between the velocity field response and this nonlinearity. By only 
considering disturbances that are periodic in the wall-parallel directions and in time, we 
restrict our attention to propagating modes. This permits us to compare our results to 
the typical statistical decomposition of point and field experimental measurements and 
simulations for both temporal and spatial data. The linearity of the equation means that 
the different frequencies and wavenumbers are only coupled through the forcing. 

Of necessity, there is a steady component that varies only in the wall-normal direction, 
which we identify as the turbulent mean profile. This yields for the perturbations a linear 



A critical layer model for turbulent pipe flow 



37 



operator identical to that obtained by linearisation about the turbulent mean profile. 
However, we do not assume small perturbations, and such a linearisation would set the 
forcing to zero. 

The mean profile is sufficient information to calculate the resolvent at any other 
wavenumber-frequcncy combination, and we use a mean profile determined from the 
Princeton superpipe experiments. 

Once the resolvent is known, we find the forcing shape at a wavenumber-frequency 
combination that gives the largest velocity field response, under the assumption that it 
is likely to be a more dominant structure in turbulent flows than other velocity field 
responses associated with smaller amplification. 

The results of the model lend themselves to an interpretation, at each wavenumber- 
frequency combination, involving propagating modes with convection or phase velocity 
equal to the velocity somewhere in the interior of the flow. Building on the analysis of 
the Orr-Sommerfeld equation in unstable flows, we identify two regions in which the 
action of viscosity is required for each wavenumber-frequency combination. They are in 
the vicinity of the critical layer, where the phase velocity is equal to the local mean 
velocity, and close to wall, where the boundary condition must be met. We designated 
these the "critical" and "wall" layers. By analogy with the upper and lower branches of 
the neutral curve in unstable flow, there exist two relative wall-normal locations for these 
layers: an upper branch type configuration in which the critical layer resides at a distance 
2/+ ~ and is distinct from the the wall layer which is centred on y+ ~ and 

a lower branch type solution in which the critical layer occurs at y+ ^ R^^^^, overlaps 
with the wall layer and has a footprint down to the wall. 

It was proposed that key features of wall turbulence can be represented in the frame- 
work of a range of propagating modes with different wavenumber-frequency combina- 
tions such that they move with different convcctive velocities. The dominant modes that 
emerge from the resolvent model provide a descriptive basis for the decomposition of 
the turbulent velocity field and capture some known features of wall turbulence. In some 
sense, the conclusions of this research represent a return to the wave-like concepts of tur- 
bulence research in the mid twentieth century, but with the interpretation that coherent 
structure is a consequence of the assembly of modes and their relative motions. 

While previous researchers have identified wall layer scaling for the Reynolds stress 
(Sreenivasan et al. , Klewicki et al. ) and the importance of propagating waves (Sirovich 
et al. ) in wall turbulence, our extended formulation suggests a picture of organised com- 
plexity in wall turbulence that appears to be consistent with a range of results from 
observations of full flows. In particular, the critical layer arguments and scaling that are 
presented here lead to inner scaling of the small scales near the wall and apparent outer 
scaling of the fluctuations in the core of the pipe, as well as a mechanism for the origin 
and scaling of the very large scale motions in pipe flow, as the critical layer solutions of 
the forced linear equation. In addition, the interpretation of the wall layer as the solution 
which meets the wall boundary condition at each wavenumber-frequency combination 
implies that there must be an interaction between the critical and wall modes, an inter- 
action that has been identified in recent literature as an amplitude modulation of the 
near-wall turbulence by the very large scales (Mathis et al. ) or the phase relationship 
between the small and large scales (Chung & McKeon). The broader implications of this 
scaling include a framework with which to interpret flow over rough walls, which can be 
considered as a modification to the boundary condition that the wall layer must meet. 
This is consistent with observations in the literature that while the near-wall flow is 
disturbed by the roughness, the VLSMs are still observed (e.g. Guala et al. , 2009). 

A further observation concerns the appropriate scaling properties at the very large 
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scales. It is clear that the VLSM wavelength scales on the outer lengthscalc, in our case the 
pipe radius. The critical layer scaling requires that the velocity scale should be the cen- 
treline velocity, Ucl, which is in agreement with arguments of Jimenez et al. concerning 
the "global" modes in channel flow. The model suggests that the dominant mode with 
streamwise and spanwise wavenumbers representative of the experimental observations 
of the VLSMs in pipe flow should travel with a phase velocity equal to two-thirds of the 
centreline velocity, such that the critical layer occurs where U{y) = 2/3Ucl- This occurs 
within the overlap region in the mean velocity for all Reynolds numbers considered here, 
and as such represents a mechanism by which the effects of viscosity can extend well be- 
yond the immediate vicinity of the wall. This is in agreement with several recent results 
including studies of the mean velocity profile and mean momentum equation. 

The lack of information in this analysis on the relative phase between modes precludes 
immediate identification of structural features such as those we would expect from ob- 
servations of the real flow, including hairpin-like vortices and the near-wall cycle, etc. 
However these may be construed, at least in part, to be a consequence of the dominant, 
energetic modes moving relative to each other with different convective velocities. This 
is a subject of ongoing investigation. 

Extrapolation of the scaling arguments to other flows suggests a fundamental difference 
between internal and external flows (or at least pipes and boundary layers). While the 
influence of spatial development in the latter case has not yet been taken into account, 
preliminary observations suggest that the upper branch type critical layer solution at 
VLSM-like a wavenumber pair may dominate, such that the VSLM energy peak scales 
with ^ R^^/^ instead of the scaling observed in pipe flow. Experimental 

confirmation of the critical layer scaling proposed here is ongoing. 

Finally, the observation that the wall layer must meet the boundary condition suggests 
that the outlook is optimistic for wall-based control of friction drag. 
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Appendix A: Interpretation of resolvent analysis 

In addition to the type of transfer function analysis presented here, our resolvent 
analysis can be understood in terms of pseudospectra analysis (Trefethen & Embree 
2005). 

The e-pseudospectrum of C, Ag is defined as 

A,(/:) = {A e C : (A/-£)-i ^ e-i} . (6.1) 

That is, the e-pseudospectrum represents the spectrum of the operator C under a 
perturbation with magnitude e. This is then interpreted as bounds on the complex plane 
within which the spectrum of the perturbed operator lies. 

The level curves of the pseudospectrum therefore correspond to the level curves of the 
resolvent and the spectrum of C represents the subset of the field of complex numbers 
where e = and the resolvent is unbounded. The norm of the resolvent (||(A/ — £)~"'^||2) 
is equal to at any particular point A in the complex plane. 
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Figure 17. Spectrum, pseudospectrum and 
turbance with (k,n) = (1,10) in laminar 
(sI-C)-^ = 10^10^10^ etc. 



resolvent norm for harmonic forcing, for a dis- 
pipe flow at Re = 5 x 10"^. Level curves at 



This is illustrated in the left pane of figure 17, in which pseudospectrum of a repre- 
sentative mode in laminar pipe flow is compared with the familiar spectrum. Where the 
pseudospectrum intrudes significantly on the right-half plane, significant perturbation 
growth with a nominally stable operator can result (Meseguer & Trefethen 2003). 

For pipe flow, the spectrum of C lies in the left- half plane (£ is always stable). Since 
the maximum value of an analytic function on a region of the complex plane lies on 
that region's boundary (the maximum modulus principle), the maximal norm of the 
resolvent taken over the right half of the complex plane of a stable operator lies on the 
the imaginary axis. This corresponds to the system response to harmonic forcing and 
can be visualised by taking a slice of a pseudospectrum contour plot along the imaginary 
axis, as shown in the right-hand panel of figure 17. 

Our resolvent analysis may be understood as an attempt to characterise system non- 
normality. Performing this characterisation by formulating the initial value problem is 
popular in fluid mechanics. The two approaches are related, but not trivially. It is a 
consequence of the Hille-Yosida theorem (Curtain & Zwart 1995) that the resolvent 
obeys 

\\{xi-c)-'h^ 



q — d 

where q is the real part of A, and M and d are real numbers and 

||e^*||2 = Me*. 



M characterises the transient behaviour, and d the asymptotic response. It is an im- 
portant open problem in control theory how to analytically determine M for a given d 
(Blondcl & Megretski 2004) . The properties of the resolvent in turbulent flow have been 
studied in the context of the initial value problem in wall shear flows (Butler & Farrell 
1992a; del Alamo & Jimenez 2006; Cossu et al. 2009) and the response of turbulent 
channel flow to stochastic forcing (Farrell & loannou 1998). 
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Figure 18. Spectrum and resolvent norm of (turbulent) jC for a disturbance with (fc, n) — (1, 10) 
at _Re — 7.5 x 10*. In the left plane, the black dots show the spectrum calculated with a resolution 
of A'^ = 60 and the red crosses the same calculation for A'^ = 200. The right pane shows the 
resolvent norm across the imaginary axis. Although the spectrum is not well resolved overall 
in either case, the spectrum near the imaginary axis is well resolved quickly and the resolvent 
norm is accurate for much lower TV than is required for a full spectral analysis. 



Appendix B: Robustness of the resolvent norm to numerical and 
modelling error 

Figure 18 shows the spectrum of the operator C for the wavenumber pair (fc, n) = 
(1, 10) in turbulent pipe flow for two different spatial resolutions, = 60 and N = 200. 
We should expect that features prominent in real flows should be robust to operator 
perturbation (such as under-resolution) and as such have a high response to forcing. 
Inspection of figure 18 demonstrates this to be the case. 

In the left plane, the black dots show the spectrum calculated with a resolution of 
A'^ = 60 and the light crosses the same calculation for A^ — 200. The right pane shows 
the resolvent norm across the imaginary axis. Although the whole spectrum is not well 
resolved overall at iV = 60, the spectrum near the imaginary axis is well resolved, sug- 
gesting that the dominant part of the spectrum is resolved with low N, and both the 
resolvent norm and the mode shape (not shown) are accurate for much lower A^ than 
is required for a full spectral analysis. However high resolution is required to observe 
near-wall, inner scaling modes at high Reynolds number. These two constraints may be 
incompatible and limit the extension of these results to higher Reynolds numbers without 
additional consideration of the numerical technique. 
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